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Proton  exchange  membrane  (PEM)  fuel  cell  performance  is  directly  related  to  the  flow  channel  design 
on  bipolar  plates.  Power  gains  can  be  found  by  varying  the  type,  size,  or  arrangement  of  channels.  The 
objective  of  this  paper  is  to  present  two  new  flow  channel  patterns:  a  leaf  design  and  a  lung  design.  These 
bio-inspired  designs  combine  the  advantages  of  the  existing  serpentine  and  interdigitated  patterns  with 
inspiration  from  patterns  found  in  nature.  Both  numerical  simulation  and  experimental  testing  have  been 
conducted  to  investigate  the  effects  of  two  new  flow  channel  patterns  on  fuel  cell  performance.  From  the 
numerical  simulation,  it  was  found  that  there  is  a  lower  pressure  drop  from  the  inlet  to  outlet  in  the 
leaf  or  lung  design  than  the  existing  serpentine  or  interdigitated  flow  patterns.  The  flow  diffusion  to  the 
gas  diffusion  layer  was  found  be  to  more  uniform  for  the  new  flow  channel  patterns.  A  25  cm2  fuel  cell 
was  assembled  and  tested  for  four  different  flow  channels:  leaf,  lung,  serpentine  and  interdigitated.  The 
polarization  curve  has  been  obtained  under  different  operating  conditions.  It  was  found  that  the  fuel  cell 
with  either  leaf  or  lung  design  performs  better  than  the  convectional  flow  channel  design  under  the  same 
operating  conditions.  Both  the  leaf  and  lung  design  show  improvements  over  previous  designs  by  up  to 
30%  in  peak  power  density. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Polymer  electrolyte  membrane  (PEM)  fuel  cells  use  pure  hydro¬ 
gen  as  a  fuel  and  either  pure  oxygen  or  air  as  an  oxidant  to  convert 
chemical  energy  of  fuel  directly  into  electrical  energy  through  an 
electrochemical  reaction.  PEM  fuel  cells  are  the  popular  option  for 
alternative  power  supply  to  the  transportation  system  owing  to 
their  high  power  density,  low  operating  temperature,  low  emission 
and  compact  design. 

The  average  electric  power  from  a  single  PEM  fuel  cell  is  lim¬ 
ited  to  around  0.5  W cm-2  [1].  Hence  several  single  PEM  fuel  cells 
must  be  stacked  to  make  a  fuel  cell  stack  in  order  to  create  the 
desired  power  for  various  applications.  To  make  a  fuel  cell  stack, 
the  components  called  bipolar  plates  are  used  to  electrically,  ther¬ 
mally,  and  physically  connect  the  fuel  cells  and  increase  the  power 
output.  Another  important  feature  of  bipolar  plates  is  that  bipo¬ 
lar  plates  contain  flow  channels  that  distribute  the  fuel  or  oxidant 
evenly  over  the  membrane  electrode  assembly  (MEA).  Meanwhile, 
the  liquid  water  produced  during  the  electrochemical  reaction  will 
be  removed  from  the  cell  through  the  flow  channels  on  the  bipolar 
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plate  [2].  Therefore,  the  flow  channel  pattern  on  the  bipolar  plate 
plays  a  very  important  role  in  maintaining  good  performance  of  fuel 
cells. 

Several  existing  flow  pattern  designs  include  the  straight  or  par¬ 
allel,  single  serpentine,  triple  serpentine,  and  interdigitated  designs 
[3].  Each  design  has  its  own  advantages  and  disadvantages  that 
affect  the  fuel  cell’s  performance  differently.  Much  research  work 
has  been  conducted  to  investigate  the  effects  of  various  flow  field 
designs  on  the  performance  of  fuel  cells  [3-29].  The  design  of  the 
flow  channel  on  bipolar  plates  mainly  focuses  on  three  aspects:  the 
flow  channel  distribution  pattern  [6-13],  the  cross-section  shape 
[14-17]  and  the  land  and  channel  dimensions  [15-22].  The  gas  is 
easily  distributed  in  the  straight  or  parallel  fields.  However,  if  a 
blockage  should  occur  in  one  of  the  channels,  the  gas  has  alternate 
flow  paths  and  therefore  will  not  build  up  the  pressure  needed  to 
force  out  the  blockage  [3,8].  The  reactant  distribution  inside  the  cell 
is  therefore  not  uniform  either,  which  will  result  in  a  non-uniform 
local  current  density  distribution.  To  solve  this  blockage  problem, 
a  single  serpentine  field  is  suggested.  As  the  name  implies,  there 
is  only  one  channel  for  the  gas  to  flow,  therefore  forcing  out  water 
droplets  and  impurities.  The  entire  channel  length,  however,  is  long 
compared  to  its  cross-sectional  area,  thus  creating  a  large  pressure 
loss  across  the  field.  This  leads  to  uneven  diffusion  into  the  gas 
diffusion  layer  (GDL)  and  can  affect  the  long  term  performance  of 
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Fig.  1.  (a)  New  bio-inspired  leaf  flow  pattern,  (b)  New  bio-inspired  lung  flow  pattern. 


the  cell.  On  the  other  hand,  the  pressure  difference  between  the 
adjacent  channels  causes  the  flow  crossover  from  one  channel  to 
another  channel  which  thus  increases  the  convection  of  reactant 
in  the  GDL  [23-27].  To  improve  the  overall  pressure  loss  in  a  sin¬ 
gle  serpentine  channel,  a  triple  serpentine  can  be  used  [11  ].  In  this 
design,  three  serpentine  paths  are  connected  at  the  inlet  and  outlet. 
However  with  the  three  channels,  the  pressure  needed  to  remove 
water  droplets  or  impurity  blockages  may  not  build  up. 

Another  popular  flow  channel  is  the  interdigitated  pattern 
which  uses  a  non-continuous  set  of  inlet  and  outlet  channels  to 
force  the  gas  through  the  GDL  [9,28-29].  In  interdigitated  designs, 
the  flow  into  the  GDL  is  forced  convection  rather  than  diffusion.  This 
helps  avoid  water  build  up  and  increases  the  limiting  current  den- 


Outlet 

Fig.  2.  Subdomains  for  interdigitated  flow  field  used  in  simulation. 


sity  and  maximum  power  by  overcoming  the  limitations  of  diffusive 
mass  transfer.  The  disadvantages,  however,  is  that  higher  inlet  pres¬ 
sures  are  required  for  operation,  which  leads  to  a  high  parasitic 
power  loss  for  the  fuel  cell  system  to  operate.  Different  from  the 
aforementioned  conventional  designs,  Morgan  Fuel  Cell  also  devel¬ 
oped  biomimetic  plates,  which  show  improvements  in  peak  power 
by  up  to  16%  using  their  biomimetic  plates  [3]. 

Much  research  is  being  done  on  flow  fields  and  there  are  still 
opportunities  to  advance  knowledge  from  its  current  state.  The  goal 
of  this  research  is  to  develop  new  designs  that  combine  the  advan¬ 
tages  of  the  existing  patterns  with  inspiration  from  natural  flow 
fields.  In  this  investigation  the  natural  flow  fields  chosen  for  two 
new  designs  were  a  lung  layout  and  a  leaf  layout:  a  combination 
of  the  serpentine  and  interdigitated  pattern  geometries  shown  in 
Fig.  la  and  b. 

The  objective  of  this  paper  is  to  evaluate  the  performance  of  two 
new  bio-inspired  flow  pattern  designs  for  bipolar  plates  in  PEM  fuel 
cells.  A  computational  fluid  dynamics  (CFD)  analysis  was  used  to 
model  the  flow  through  the  newly  designed  channels  and  diffu¬ 
sion  into  the  GDL.  The  resulting  pressure  distribution  and  pressure 
loss  is  shown  visually  on  the  exchange  boundary  between  the  flow 
channel  and  GDL.  Experimentally,  the  performance  of  the  two  exist¬ 
ing  designs  and  the  new  lung  and  leaf  patterns  were  determined 
using  a  Fuel  Cell  Technologies  test  station.  These  bio-inspired  plates 
demonstrate  flow  improvements  as  well  as  improved  cell  perfor¬ 
mance.  The  new  designs  show  improved  pressure  distributions, 
decreased  pressure  loss,  and  increased  cell  performance  in  current 
density  and  peak  power  density  compared  to  the  single  serpentine 
and  interdigitated  designs. 


2.  Computational  model 

2.1  Model  assumptions  and  governing  equations 


Each  design  to  be  studied,  including  the  single  serpentine,  inter¬ 
digitated,  lung,  and  leaf  flow  patterns,  was  physically  modeled  in 
Unigraphics  NX4  and  then  imported  as  a  parasolid  to  COMSOL 
Multiphysics  3.4.  The  flow  is  assumed  to  be  steady,  laminar,  and 
incompressible  in  both  flow  channels  and  the  gas  diffusion  layer. 
The  fuel  cell  is  assumed  to  operate  isothermally.  The  momentum 
transport  in  a  flow  channel  is  governed  by  Navier-Stokes  Eq.  (1), 
and  the  momentum  transport  in  the  porous  gas  diffusion  layer  is 
governed  by  Brinkman’s  Eq.  (2).  The  mass  transport  in  both  regions 
is  governed  by  the  continuity  Eq.  (3). 

puxVu  =  V[-PI  +  ft(Vu  +  VuT)]  (1) 


( v/k)u  =  V 


-PI  + 


u(Vu  +  VuT) 


(2) 


V  x  u  =  0 


0) 
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where  p  is  the  fluid  density,  kgrrr3,  u  is  the  velocity,  ms-1,  P 
is  pressure,  Pa,  /z  is  the  fluid  kinematic  viscosity,  ms-2,  /c  is  the 
permeability  of  gas  diffusion  layer,  m2,  and  £  is  the  porosity  of  the 
gas  diffusion  layer. 

2.2.  Model  geometry  and  boundaries 

In  both  the  flow  field  and  GDL  subdomains,  air  is  used  as  the 
working  fluid.  A  constant  temperature  of  75  °C  and  the  calculated 
pressure  were  used  to  establish  the  density  and  dynamic  viscosity 
of  the  fluid.  The  GDL  subdomain  was  connected  to  the  flow  field 
through  a  continuity  boundary  on  the  exchange  surface  between 
the  GDL  and  flow  field.  For  Brinkman’s  Eq.  (2),  a  porosity  of  0.78 
and  permeability  of  lE-8m2  were  used.  Barbir  has  reported  a 
permeability  of  5-10  x  10-12  m2  for  carbon  paper  GDLs  used  in 
other  simulations  [30],  however,  the  permeability  in  this  model 
was  increased  to  enhance  the  flow  visualization  though  the  GDL 
exchange  surface.  A  GDL  thickness  of  1.016  mm  was  also  used  to 
enhance  the  flow  visualization  through  the  GDL  exchange  surface. 
The  three  subdomains  (inlet  flow  field,  outlet  flow  field,  and  GDL) 
are  shown  in  Fig.  2  color  coded  by  the  legend  for  an  interdigitated 


Table  1 

Input  variables  for  computational  model. 


Variable  name 

Value 

Description 

T(°C) 

75 

Temperature 

Po (Pa) 

0 

Initial  pressure 

E 

0.78 

Porosity  of  GDL 

k  (m2) 

IE— 8 

Permeability  of  GDL 

fgdi  (mm) 

1.016 

Thickness  of  GDL 

Vi  (ms-1) 

20.75 

Inlet  velocity 

Po (Pa) 

0 

Outlet  pressure  (gage  pressure) 

wc  (mm) 

0.7874 

Flow  field  channel  width 

dc  (mm) 

1.016 

Flow  field  channel  depth 

wr  (mm) 

0.8 

Rib  width 

pattern  design.  A  similar  arrangement  was  used  for  the  single  ser¬ 
pentine,  leaf,  and  lung  designs. 

Four  main  boundary  conditions  were  created  for  each  flow  pat¬ 
tern:  inlet,  outlet,  walls,  and  GDL  exchange.  In  the  inlet  boundary, 
the  fluid  velocity  was  given.  In  the  outlet  boundary,  the  outlet  pres¬ 
sure  was  set  to  atmospheric  pressure.  A  ‘no-slip’  condition  was  used 
for  all  of  the  wall  boundaries.  Finally,  the  Navier-Stokes  Eq.  (1)  for 
the  flow  in  the  channel  and  Brinkman’s  Eq.  (2)  for  the  flow  in  porous 
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Fig.  3.  (a)  Pressure  distribution  in  a  single  serpentine  pattern,  (b)  Pressure  distribution  in  an  interdigitated  flow  pattern,  (c)  Pressure  distribution  in  a  leaf  flow  pattern,  (d) 
Pressure  distribution  in  a  lung  flow  pattern. 
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Fig.  4.  Pressure  distribution  and  velocity  profile:  the  cut  view  in  the  middle  of  each  design  pattern. 


gas  diffusion  layer  are  coupled  along  the  interface  between  the  GDL 
and  channel  through  a  ‘continuity  boundary’  condition  given  as 
bellows: 

^NS  —  ^Brinkman  ^NS  —  ^Brinkman  (4) 

Table  1  summarizes  values  and  constants  used  in  the  computa¬ 
tional  model. 


2.3.  Simulation  results  and  discussion 

The  SPOOLES  solver  (convergence  tolerance  =  IE-3  and  maxi¬ 
mum  iterations  =  100)  in  Comsol  Multiphysics  3.4  was  used  to  solve 
these  equations  for  pressure  and  velocities  in  the  three  Cartesian 
axis  directions.  Pressure  loss  from  the  GDL  input  to  output  and  the 
velocity  normal  to  the  GDL  exchange  surface  were  computed  to 
assess  the  four  designs. 
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Table  2 

The  computed  pressure  loss  for  various  flow  patterns. 

Flow  field 

Pressure  loss  (Pa) 

Single  serpentine 

41.23 

Interdigitated 

38.023 

Leaf 

25.794 

Lung 

27.325 

The  actual  pressure  drop  the  GDL  input  at  the  channel  inlet  to 
the  output  at  the  channel  outlet  is  shown  in  Table  2.  Fig.  3  shows  the 
pressure  distribution  in  different  flow  patterns.  Both  the  leaf  and 
lung  patterns  had  an  improved  pressure  loss  over  the  single  ser¬ 
pentine  and  interdigitated  patterns.  The  values  of  the  pressure  drop 
for  the  single  serpentine  and  interdigitated  patterns  are  consistent 
with  simulations  performed  by  Karvonen  et  al.  [31].  The  pressure 
distributions  along  the  GDL  exchange  were  also  more  uniform  in  the 
new  designs  as  illustrated  in  Fig.  3c  and  d.  Both  the  leaf  and  lung 
designs  result  in  less  high  pressure  regions  (shown  in  red)  which 
help  create  a  uniform  diffusion  into  the  GDL.  The  pressure  distribu¬ 
tion  and  the  corresponding  velocity  profile  were  further  examined 
by  a  cut  view  in  the  middle  of  each  flow  pattern  design  as  shown 
in  Fig.  4.  A  large  magnitude  of  velocity  is  observed  in  both  interdig¬ 
itated  and  serpentine  flow  channels  due  to  a  large  pressure  drop. 


The  velocity  in  both  leaf  and  lung  designs  are  relatively  small  but 
more  uniform.  Also,  the  pressure  drop  in  both  leaf  and  lung  designs 
in  the  flow  direction  is  much  smaller  than  that  in  interdigitated  and 
serpentine  designs.  The  Z  velocity  contour  profile  in  Fig.  5a-d,  show 
that  the  diffusion  through  the  GDL  exchange  surface  is  more  uni¬ 
form  in  the  new  leaf  and  lung  designs  and  there  are  minimal  spikes 
or  dead  zones  that  negatively  affect  the  cell’s  performance.  A  high 
and  uniform  diffusion  flow  rate  to  the  GDL  will  improve  the  fuel 
cell  performance,  which  can  be  seen  in  next  section. 

3.  Experimentation 

In  order  to  check  the  performance  of  new  designs,  four  single 
PEM  fuel  cells  are  assembled  with  the  serpentine  flow  pattern, 
interdigitated  flow  pattern,  leaf  flow  pattern  and  lung  flow  pattern 
on  the  bipolar  plates,  respectively.  The  Nil 2  membrane  electrode 
assembly  (MEA)  from  FuelCellStore.com  was  used  for  all  the  four 
designs  with  an  anode  and  cathode  loading  of  0.5  mgPtcm-2  each. 
This  MEA  performed  well  and  produced  consistent  results  over  the 
extensive  weeks  of  testing.  In  order  to  fabricate  bipolar  plate  chan¬ 
nels,  a  Computer-Controlled  Bridgeport  mill  at  Oakland  University 
was  used  to  cut  the  patterns  on  the  high  purity  POCO®  electro¬ 
graphite  plates.  In  order  to  avoid  any  MEA  damage  or  deformation, 
the  bolt  torques  for  each  assembly  were  kept  at  a  constant  120  in.  lb 


x  (m)  x  (m) 

Fig.  5.  (a)  Contour  profile  of  the  Z  velocity:  interdigitated  flow  pattern,  (b)  Contour  profile  of  the  Z  velocity:  serpentine  flow  pattern,  (c)  Contour  profile  of  the  Z  velocity:  leaf 
flow  pattern,  (d)  Contour  profile  of  the  Z  velocity:  lung  flow  pattern. 
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Fig.  6.  (a)  Leaf  design:  effects  of  operation  temperature-polarization  curve,  (b)  Leaf  design:  effects  of  operation  temperature-power  density,  (c)  Lung  design:  effects  of 
operation  temperature-polarization  curve,  (d)  Lung  design:  effects  of  operating  temperature-power  density. 


Fig.  7.  (a)  Leaf  design:  effects  of  relative  humidity-polarization  curve,  (b)  Leaf  design:  effects  of  relative  humidity-power  density,  (c)  Lung  design:  effects  of  relative 
humidity-polarization  curve,  (d)  Lung  design:  effects  of  relative  humidity-power  density  curve. 
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Table  3 

Summary  of  testing  conditions.  For  all  test  groups:  anode:  99.9999%  pure  hydrogen; 
300  seem;  cathode:  compressed  air:  1000  seem. 


Temperatures 

(°C) 

Back  pressure 
(atm) 

Relative  humidity 
(anode:cathode) 

Test  Group  1 :  (temperature)  35 

1 

100%:  100% 

55 

1 

100%:100% 

75 

1 

100%:100% 

85 

1 

100%:100% 

Test  Group  2:  (pressure)  75 

0 

100%:  100% 

75 

1 

100%:100% 

75 

2 

100%:100% 

Test  Group  3:  (RH)  75 

1 

100%:0% 

75 

1 

100%:30% 

75 

1 

100%:60% 

75 

1 

100%:100% 

using  a  beam  deflection  torque  wrench.  Previous  tests  [32]  found 
this  bolt  torque  value  to  be  safe  and  high  performing. 

The  University  Model  Test  Station  by  the  Fuel  Cell  Technologies 
Inc.  was  used  to  accurately  measure  the  current  density  and  voltage 
during  testing.  This  test  station  maintains  the  inlet  reactant  tem¬ 
perature,  cell  temperatures,  outlet  backpressures,  and  inlet  relative 
humidity.  Except  for  the  pressure  regulators  at  the  inlet,  all  of  these 
conditions  are  altered  and  monitored  using  LabVIEW  and  recorded 
to  a  data  sheet  output. 

Three  test  groups  were  performed  to  find  optimum  conditions 
for  the  new  lung  and  leaf  designs:  cell  temperature,  relative  humid¬ 
ity,  and  backpressures.  In  all  of  the  groups,  the  anode  and  cathode 
pre-heaters  were  set  to  5°C  higher  than  the  anode  and  cath¬ 
ode  humidifier  temperatures  to  avoid  condensation  before  the  gas 
enters  the  fuel  cell.  During  the  backpressure  tests,  the  tank  regu¬ 
lators  were  increased  from  30psi  to  60psi  in  order  to  provide  a 
sufficient  pressure  difference  between  the  inlet  and  outlet.  The  test 
groups  were  also  run  in  successive  order  with  a  cell  purge  between 
each  test.  For  example  during  the  temperature  tests,  the  first  test 
was  at  the  lowest  temperature  while  the  final  test  was  at  the  high¬ 


est  temperature.  A  similar  progression  was  used  for  the  relative 
humidity  (0-100%)  and  backpressure  (0-2  atm)  experiments.  The 
operating  conditions  summarized  in  Table  3,  were  repeated  for  each 
flow  field  design.  The  constant  flow  rates  of  300  seem  ultra  pure 
hydrogen  and  1000  seem  compressed  air  are  supplied  to  the  anode 
and  cathode  electrode,  respectively.  A  polarization  plot  of  voltage 
(V)  vs.  current  density  (A cm-2)  and  a  power  density  (Wcnrr2)  vs. 
current  density  (A cm-2)  plot  for  all  four  patterns  were  used  to 
show  which  design  achieves  the  best  performance  and  peak  power 
density. 

4.  Experimentation  results 

4.1.  Temperature  effects 

Fig.  6  shows  the  effects  of  operation  temperature  on  the  fuel 
cell  with  the  new  leaf  and  lung  pattern  design.  The  new  leaf  and 
lung  designs  show  improvements  in  performance  as  temperature 
increases.  This  performance  improvement  peaks  at  65-75  °C  and 
then  declines  during  the  85  °C  test.  One  explanation  for  this  is  that 
at  temperatures  greater  than  75  °C,  the  incoming  water  vapor  does 
not  hydrate  the  MEA  well,  but  rather  dries  it  out  which  decreases 
the  performance.  The  reason  is  that  the  thermal  couple  in  our  test 
rig  was  located  in  the  cathode  end  plate.  Therefore,  the  real  tem¬ 
perature  in  the  membrane  of  the  fuel  cell  should  be  a  few  degrees 
higher  than  the  actual  reading  from  the  thermal  couple  due  to  heat 
loss  from  the  metal  end  plates.  During  experimentation,  it  was  also 
noticed  that  the  gas  at  the  exit  was  mostly  water  vapor  at  higher 
temperatures  and  more  visible  water  droplets  at  lower  tempera¬ 
tures.  Since  less  water  remains  in  the  GDL  at  higher  temperatures, 
the  MEA  can  become  dehydrated  which  thus  decreases  the  perfor¬ 
mance  of  cell. 

4.2.  Relative  humidity  effects 

Fig.  7  shows  the  effects  of  relative  humidity  at  the  inlet  of  the 
fuel  cell  on  the  cathode  side.  The  new  lung  and  leaf  pattern  designs 


Fig.  8.  (a)  Leaf  design:  effects  of  operating  pressure-polarization  curve,  (b)  Leaf  design:  effects  of  operating  pressure-power  density,  (c)  Lung  design:  effects  of  operating 
pressure-polarization  curve,  (d)  Lung  design:  effects  of  operating  pressure-power  density. 
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Fig.  9.  (a)  Polarization  curve  comparison  for  various  flow  pattern  designs:  P=  2  atm, 
T=  75  °C,  RH  =  100%.  (b)  Power  density  comparison  for  various  flow  pattern  designs: 
P= 2  atm,  T=  75  °C,  RH  =  100%. 


follow  a  similar  trend  with  an  increase  in  cathode  relative  humid¬ 
ity.  With  a  higher  humidified  oxidant,  the  membrane  stays  well 
hydrated,  which  facilitates  the  movement  of  electron  from  anode 
to  cathode  electrode.  As  a  result,  there  is  a  low  ohmic  loss  in  the 
membrane  as  the  humidity  increases.  A  high  current  density  is  also 
achieved  as  the  humidity  increases.  However,  this  trend  also  has 
a  saturation  point  in  the  regions  where  water  is  readily  produced 
from  the  high  current  demand.  Water  build  up  can  occur;  how¬ 
ever,  because  it  reaches  the  highest  peak  power  density  during  the 
100%  humidity  condition,  this  remains  the  optimum  condition.  This 
trend  can  be  observed  in  Fig.  7a  and  b  for  the  Leaf  and  Lung  designs. 
The  lung  design  shows  a  better  performance  than  the  leaf  design 
under  the  same  relative  humidity.  The  maximum  current  density 
reached  by  the  lung  design  a  little  higher.  The  lung  design  removes 
water  more  easily  than  the  leaf  design  considering  its  symmetrical 
parallel  design  feature. 

4.3.  Backpressure  effects 

Fig.  8  shows  the  effects  of  backpressure  at  the  outlet  of  the  fuel 
cell.  During  the  test,  the  backpressure  was  maintained  to  be  the 
same  on  both  anode  and  cathode  sides.  As  pressure  increases,  the 
fuel  cell  performs  better  with  either  a  leaf  or  lung  design.  This 
increase  in  performance  is  very  obvious  when  the  backpressure 
increases  from  0  atm  to  1  atm,  and  it  is  not  so  obvious  when  the 
backpressure  increases  from  1  atm  to  2  atm.  This  is  consistent  with 
the  founding  for  the  fuel  cell  with  an  interdigitated  channel  design 
[28].  The  reaction  rate  is  increased  because  that  more  oxidant  is 
available  at  higher  pressures  and  forced  into  the  GDL  more  readily. 
However,  one  drawback  with  increasing  the  backpressure  is  that 
water  can  build  up  due  to  a  lower  pressure  differential  between  the 
inlet  and  outlet,  which  causes  unsteady  test  results.  An  optimum 
backpressure  of  2  atm  is  shown  in  Fig.  8a  and  b  for  the  leaf  and  lung 


designs;  however,  frequent  purges  were  needed  to  remove  water 
for  the  higher  reaction  rate. 

4.4.  Performance  comparison  between  the  serpentine, 
interdigitated,  leaf  and  lung  designs 

Serpentine,  interdigitated,  leaf,  and  lung  designs  are  compared 
at  the  same  operating  conditions:  backpressure  of  2  atm,  temper¬ 
ature  of  75  °C,  and  a  relative  humidity  of  100%.  As  shown  in  Fig.  9, 
the  new  leaf  and  lung  designs  show  improvement  in  overall  perfor¬ 
mance  and  peak  power  density  by  up  to  30%  compared  to  either  the 
serpentine  or  interdigitated  design.  The  leaf  design  has  the  best  per¬ 
formance  followed  by  the  lung,  interdigitated  and  single  serpentine 
design. 

5.  Conclusions 

In  PEM  fuel  cells,  bipolar  plates  are  major  components  that 
physically,  electrically,  and  thermally  connect  the  cell  together  to 
produce  electrical  power.  Flow  patterns  on  these  bipolar  plates  sig¬ 
nificantly  affect  the  immediate  and  long  term  performance  of  the 
fuel  cell.  This  research  shows  promising  results  for  bio-inspired  flow 
patterns  with  noticeable  improvements  in  pressure  loss,  overall 
performance,  and  peak  power  density.  The  optimum  operating  con¬ 
ditions  of  65-75  °C  cell  temperature,  2  atm  backpressure,  and  100% 
relative  humidity  are  found  for  the  new  leaf  and  lung  designs.  Under 
these  operating  conditions,  the  leaf  design  has  the  best  perfor¬ 
mance  of  the  two.  Both  the  leaf  and  lung  design  show  improvements 
over  previous  designs  by  up  to  30%  in  peak  power  density. 

To  achieve  further  improvements,  the  following  should  be  inves¬ 
tigated:  tapered  ends,  rounded  corners,  and  surface  treatment  of 
the  plates.  Tapered  ends  may  further  improve  the  pressure  loss  on 
each  channel  as  shown  in  previous  research.  Although  milling  may 
be  difficult,  rounded  corners  may  help  create  more  uniform  flow 
through  the  channels  and  thus  more  consistent  diffusion  into  the 
GDL.  Heat  surface  treatment  of  the  plates  will  create  a  more  uni¬ 
form  surface  and  possibly  lower  the  contact  resistance  between 
the  MEA  and  bipolar  plate.  Future  work  with  these  PEM  fuel  cells 
includes  a  transparent  endplate  fuel  cell  in  order  to  visualize  the 
water  build  up  inside  the  cell.  This  transparent  fuel  cell  will  clearly 
show  any  areas  of  water  or  impurity  build  up  and  suggest  possible 
improvements  to  the  designs  tested. 
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